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Abstract 

In detection of electrons from /i — > e conversion process the 
monochromatic electron spectrum is transformed due to a photon emis- 
sion and fluctuations of energy loss in a target. The selection criterion 
of /Lt — > e{"f) conversion events is an electron momentum above the 
threshold momentum of 103.5 MeV/c , which corresponds to the max- 
imum energy loss of 1.5 MeV. Radiative corrections including a virtual 
photon correction and soft photon emission below 1.5 MeV lead to a 
reduction by about 10% in the probability of /i — > e conversion process 
calculated without radiative corrections. 

The soft photons emission below 1.5 MeV contributes to a change of 
electron spectrum from monoenergetic one at 105 MeV to a spectrum 
with a low energy tail for the process of /i — > 6(7) conversion. However 
the effect of smearing of the initial momentum distribution due to 
the soft photon emission is small in comparison with a smearing due 
to energy loss fluctuations in a target. The average energy of soft 
photons emitted below 1.5 MeV is found to be 40 keV. The soft photon 
approximation is a good description for /i — > 6(7) conversion process 
with photons emitted below 1.5 MeV. 



*Permanent address: Institute for Nuclear Research, 60-th Oct. pr. 7a, Moscow 
117312, Russia 
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Introduction 



In this work we study radiative corrections to the process oi ii —> e con- 
version, i.e. we consider the process oi ji ^ 6(7) conversion including the 
elastic process of /j, ^ e conversion and the process of — > e conversion with 
a photon emission below a certain energy ujmax- The fi ^ e conversion is 
the process of a transformation of a muon into an electron in the presence 
of a nucleus. The clear signature of /x — > e conversion is an appearance in 
the final state of a monochromatic electron of the energy 

Emax = E^-^^, E^ = m^-{l-{a-Zf/2) (1) 

where Ma is a nuclear mass, m^u is a muon mass, a is the fine structure 
constant, Z is a nuclear charge. In particular, Emax is 104.963 MeV for 
aluminum. 

In an experimental setup a monochromatic spectrum of — > e conversion 
electrons at 105 MeV is transformed to a distribution with a low energy 
tail. This transformation in the electron spectrum can be either due to 
fluctuations of energy loss of electrons in a target or because of a photon 
emission in the process of jj, ^ 6(7) conversion. The selection criterion of 
/J, — > 6(7) conversion events is an electron momentum above the threshold 
momentum of 103.5 MeV/c , which corresponds to the maximum energy 
loss of 1.5 MeV. 

It is important to note that in real experiments a pure elastic process is 
not physically observable since due to the presence of charged particles in 
the initial, intermediate or final states the process can be accompanied by an 
emission of a photon. Therefore the measured probability of the process of 
— > 6(7) conversion is a sum of probabilities of the pure elastic process cor- 
rected by virtual photon exchange and of a real soft photon emission below 
a maximum energy uJmaxi which in our case is 1.5 MeV. In fact if the emitted 
real soft photon is not detected, one can not distinguish a bremsstrahlung 
process from the elastic one. 

Because the energy release in the process of /x — > 6(7) conversion is 
about AE = 105 MeV one could expect that the virtual photon corrections 
and the soft photon emission can be enhanced by a large logarithmic factor 
ln{mn/me) where nie is an electron mass. Radiative corrections lead to a 
change in the probability of |U — > e conversion process calculated without 
radiative corrections. 

The real soft photons emission contributes to a change of electron spec- 
trum from monoenergetic one at 105 MeV to a spectrum with a low energy 
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tail. However this effect of smearing of the initial momentum distribution 
due to the soft photon emission should be compared with a smearing due to 
energy loss fluctuations in a target. 

The detailed analysis of hard photon emission requires knowledge of 
specific particle physics models. In order to avoid complications in analy- 
sis of high energy photon spectrum and to understand qualitatively a hard 
photon emission the radiative pion decay is considered. By analyzing this 
decay which possesses the same kinematic properties and a comparable en- 
ergy scale as /i — > e conversion the range of application of the soft photon 
approximation can be found. 

It is worth to note that the hard photon emission with photons produced 
at energy above 1.5 MeV is not important for simulation of /Li — 6(7) con- 
version process. Hard photons carry away a significant part of muon energy 
thus reducing an electron energy but due to the selection criterion 103.5 
MeV/c on the threshold momentum such electrons are rejected and such 
events do not contribute to useful events. 

Virtual photon corrections 

Let us consider /i — > e conversion process in which a virtual photon emission 
is involved. The Feynman diagrams describing virtual photon emission [Ij 
are presented in Figure [H 




Figure 1: Virtual photon corrections to lepton vertex and ingoing and 
outgoing leptons in the process of muon-electron conversion. 

It can be shown that an integral describing the radiative correction to 
the lepton vertex is infrared divergent. This integral can be artificially well- 
defined with an introduction of an infinite small mass A of a photon. With 
this modification the leading term in leading logarithm approximation is 
given by 
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(2) 



where T^i is the — > e conversion rate calculated without radiative 
corrections. 

Soft photon emission 

Feynman diagrams describing the process of inner bremsstrahlung are pre- 
sented in Figure [21 



Figure 2: Inner bremsstrahlung in the process of muon-electron conversion. 

Let's consider the case of emission of low energy photons. In this limit we 
assume that a photon energy is much less than energies of particles taking 
part in muon-electron conversion and momentum transfered to a nucleus. 
On the other hand we assume the validity of the perturbation theory, i.e. the 
smallness of a parameter a ■ ln(m^/me) « 1 which provides a suppression 
of multi-photon emission. 

Following a standard technique (see e.g. [I] ) in the limit of photon 
low energies the photon emission can be taken into account by adding an 
external photon line to a diagram of the elastic process. In this case the most 
important are diagrams with an emission from external charged particles 
lines, in all other parts of diagrams one can neglect changes in momenta due 
to the soft photon emission. Therefore Dirac spinors of initial muon and 
outgoing electron are modified in the following way 
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where k, ki, q are the four- vectors of muon, electron and photon, re- 
spectively, £ is the polarization of a photon. A photon emission from a 
non-relativistic nucleus is neglected. 

By using Dirac equations these relations can be simplified 

(fee*) 

-e- -jj^Uf, , 

- (fcie*) 
Ue^ e-Ue - -jj^ . 

As a result amplitude M of muon-electron conversion with the soft pho- 
ton emission can be factorized according to 



M = Mei-e 



(kie*) (ke*) 



(3) 



(hq) (kq) 

where M^i is the amplitude of the same process of ^ e conversion 
without photon emission. This is well-know general result which is process 
independent. 

Summing up over the photon polarization and by applying standard 
technique of calculation of process probability the differential probability 
to emit a single soft photon with momentum q in the process of /Li — > e 
conversion is given by 
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Integration of Eq. Q over an electron momentum within kinematic limits 
gives the spectrum of soft photons in the process of muon-electron conversion 
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Note that the same expression for the spectrum of soft photon emission 
can be obtain by applying the classical consideration [2j instead of Feynman 
diagram technique. 

The limit a; ^ in Eq.Q represents a logarithmic divergence meaning 
that the total probability of emitting a single very soft photon by a charge 
particle is infinite, this is well-known infrared catastrophe. By introducing 
a small photon mass A the integral over uj becomes well-defined and the 
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probability of the soft photon emission is derived from Eq. ([5]) by integrating 
over the photon energy from the small photon mass to the maximum energy 
^max of emitted soft photons 
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Probability of radiative muon - electron conversion 

The inner bremsstrahlung spectrum diverges for low energy photons but 
this divergence in the total rate is canceled by virtual photon corrections. 
Therefore physically observable is the sum of probabilities of the elastic 
process corrected by the virtual photon emission and of the soft photon 
emission. 

By adding the probabilities Tyirt and Tgoft the infrared divergences from 
the soft bremsstrahlung and from the lepton vertex corrections cancel each 
other producing a finite probability of radiative /i — > 6(7) conversion with 
the emission of soft photon with the energy below cOmax- 
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Note that this is the probability of radiative /i 6(7) conversion calcu- 
lated in the first order in electromagnetic constant a. 

Table [1] shows the ratio of probabilities r/Tg; of /i ^ e conversion 
calculated with and without radiative corrections. Radiative corrections 
include the virtual photon emission and soft photon emission with energy 

below UJmax- 



tOmax (MeV) 


0.1 
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1.5 


2 


2.5 




0.828 


0.867 


0.885 


0.895 


0.902 


0.907 



Table 1: The ratio of probabilities of // — > e conversion calculated with and 
without radiative corrections versus maximum soft photon energy uJmax ■ 



One can see that radiative corrections reduce the probability calculated 
without radiative corrections. For ujmax = l-5MeV this corrections is about 
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11%. 



Spectrum of conversion electrons 

In this section we consider the process of soft photon emission with energy 
below iOmax = l-5MeV because as it was discussed above the selection 
criterion of /x — > 6(7) conversion events is an electron momentum above the 
threshold momentum of 103.5 MeV/c , which corresponds to the maximum 
energy loss of 1.5 MeV. 

It is important to note that one does not need to simulate the photon 
spectrum down to very small photon energies. A straggling in a target leads 
to a smearing of the electron momentum. To get the electron momentum 
distribution from a target a pattern recognition and track reconstruction 
procedure p] based on the Kalman filter technique [1] was applied taking 
into account backgrounds, delta-rays and straw inefficiency. It will be shown 
below that an electron momentum distribution from a target can be fitted 
by a Gaussian with a standard deviation a = 200 keV/c. Therefore a pho- 
ton emission below 200 keV is indistinguishable from the elastic process of 
muon - electron conversion and it is not required to simulate this photon 
emission. In the following we will simulate photon emission in the range 

from UJmin — 

WOkeV to 

^max — l.bAdeV . The particular choice of uJrnin 
is not important because as we mentioned above for photon energies below 
200 keV the process of ^ ^ 6(7) conversion will look as the elastic process. 
Moreover a formula for the total probability of — > 6(7) process contains 
logarithm of comin and therefore it is not very sensitive to oJmin- 

The relative probability Tgoft/^ei of photon emission in the range from 
'^min to oJrriax with respcct to the probability of the corresponding elastic 
process can be calculated by integrating the photon spectrum Eq.([5]) of 
radiative — > e conversion over photon energy: 



Tel 

Figure [3] shows the relative probability to emit a photon for the radiative 
/i — > e conversion. 

Table [2] shows the relative probability to emit a photon in the range from 
i^min to LOmax for radiative muon-electron conversion calculated by Eq.Q. 

Let's consider two subranges of the photon energy: below 100 keV and 
from 100 keV to 1.5 MeV. Below 100 keV we consider the process as the 
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Relative Probability of Radiative |i -> e Conversion 
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Figure 3: Relative probability to emit a photon in the range from uj^nin 
to 

^max for radiative muon-electron conversion in the approximation of soft 
photon emission. Minimal photon energy is ojmin = 100 keV . 
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probability (%) 
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5.7 


6.7 


7.4 


8.4 


9.1 


9.7 



Table 2: Relative probability to emit a photon in the range from 100 keV 

to OJmax ■ 



elastic one occurring according to Xable [T] for uJmax — 100 keV with the 
relative rate F/Fei = 82.8%. 

Photons emitted in the subrange 100 keV - 1.5 Mev carry away a part 
of energy released in ^ 6(7) conversion. This leads to the appearance 
of low energy tail in the electron distribution. According to Table [2] with 
the relative rate Tgoft/F(,i = 6.7% in the process of — > 6(7) conversion the 
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electron energy is distributed in the range 103.5 - 104.863 MeV. 

In Monte Carlo simulations the elastic process is simulated with the 
probability 82.8/(82.8 + 6.7)% « 92.5% and the process of soft photon 
emission in the range from 100 keV to 1.5 MeV with the probability 
6.7/(82.8 + 6.7)% ?a 7.5%. The total number of useful fi — > 6(7) con- 
version events including the elastic process and one photon emission is 

= (0.828 + 0.067)iVo = 0.895iVo where Nq is the number of elastic events 
calculated without radiative corrections. We take into account the photon 
energy distribution in the simulation by sampling the photon spectrum ([5]) 
according to the equation uj = ujmin • C^,'"""' Y where r is distributed uniformly 
from to 1. 

The spectrum of electrons produced in the process of radiative ^ e 
conversion is shown in Figure H] in linear (left) and logarithmic (right) 
scale. Due to the photon emission a low energy tail in the electron spectrum 
appears. 
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Figure 4: Momentum distribution of electrons produced in the process of 
/i — > 6(7) conversion in linear and logarithmic scale. Soft photon emission 
is taking into account. 

The momentum distribution of electrons from a target smeared by strag- 
gling is shown in Figure [5] in linear (left) and logarithmic (right) scale by 
neglecting the photon emission. 

By comparing of Figures H] and [5] one can see that the effect of the 
photon emission is small in comparison with the effects of straggling in the 
target. 
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Figure 5: Momentum distribution of electrons from a target. This distri- 
bution does not include a photon emission. 



In the soft photon approximation the average energy of photons emitted 
in the range from Wmin to LOr, 



-'max 



u) = L I oj — -^dui = 7^ . 9) 



me 



For uJmin ^^max = 1-5 MeV this equation gives Co ^ 0.04 MeV. 

The classical formula from [2] also describes correctly the soft photon 
emission but the extension of this formula to high energy photons is a rough 
approximation. This is especially important in calculations of the average 
radiated energy because this quantity depends crucially on the high energy 
behavior of the spectrum. In order to describe the high energy end of the 
photon spectrum specific particle physics models are required. 



Inner bremsstrahlung in the radiative pion decay 

In order to understand qualitatively a hard photon emission and to validate 
the range of application of soft photon approximation the radiative pion 
decay is considered. In previous section the spectrum of soft photons was 
found for the process of fi ^ e{j) conversion. In order to get predictions 
for the high energy part of the photon spectrum one has to relay on a vari- 
ety of specific models for the fi — e vertex. In order to avoid complications 
in analysis of the high energy photon spectrum we just consider the inner 



10 



bremsstrahlung in the process of radiative pion decay that possesses similar 
kinematic behavior and similar mass scale as /i — > 6(7) conversion. Equa- 
tions describing the photon and electron spectra for the radiative pion decay 
are given in Appendix. 

Figure O shows the differential photon spectrum in relative photon en- 
ergy X = 2a;/m7r for the inner bremsstrahlung in the radiative pion de- 
cay. Lower curve is the result of calculations based on the exact formula 
of Appendix. Upper curve is obtained in the approximation of soft photon 
emission where only a term proportional to x~^/n(m,r/™e) contributes. 
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Figure 6: Differential spectrum in relative photon energy x = 2a;/m,r for 
the inner bremsstrahlung in the radiative pion decay. Lower curve is based 
on the exact formula, upper curve is obtained in the approximation of soft 
photon emission. 

It follows from this plot that up to the photon energy w = 5 MeV the 
exact photon spectrum can be approximated by soft photon emission 
with the precision better than 10%. At high energies the spectra differ 
significantly because as x tends to the kinematic boundary x = 1 — 
where A = me/rriT^ the exact spectrum tends to 0. 

This behavior of the photon spectrum is of general nature. In particular 
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in the case of /i — > 6(7) conversion low energy spectrum is exactly the 
same if one uses the proper mass scale. At high energies one could expect 
qualitatively similar photon spectra for /x 6(7) conversion and radiative 
pion decay. 

Figure [7] shows the differential electron spectrum in relative electron 
energy y = 2i?/m,r for the inner bremsstrahlung in the radiative pion decay. 
Lower curve is based on the exact formula of Appendix, upper curve is 
obtained in the approximation of soft photon emission. This spectrum is 
divergent as y tends to 1 + because this limit corresponds to the emission 
of low energy photons. 

Electron Spectrum of 7t ^ e v y | 
■a. 
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Figure 7: Differential spectrum in relative electron energy y = 2E/m.,j for 
the inner bremsstrahlung in the radiative pion decay. Lower curve is based 
on the exact formula, upper curve is obtained in the approximation of soft 
photon emission. 

By integrating the photon spectrum of radiative pion decay over photon 
energy one gets a relative probability F/Fg; of photon emission in the range 
from oJmin to oJmax with respcct to the probability of the corresponding 
elastic process (see Appendix). Figure [8] shows the relative probability of 
photon emission for the radiative pion decay. Minimal photon energy is 
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^min = 100 keV. Lower curve is based on the exact formula, upper curve is 
obtained in the approximation of soft photon emission. 
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Figure 8: Relative probability to emit a photon in the range from uj'jjiin 
to uJmax for the radiative pion decay. Minimal photon energy is uJmin = 
100 keV. Lower curve is based on the exact formula, upper curve is obtained 
in the approximation of soft photon emission. 

In particular, it follows from Figure [8] that the relative probability to 
produce a photon in the energy range 100 keV - 1.5 MeV is 7% which is 
comparable, due to an increase in a rate of photon emission at low ener- 
gies, with the probability of 7.8% to produce a high energy photon in much 
broader range 1.5 MeV - 70 MeV. Note that the range 100 keV - 1.5 MeV is 
described with high precision (better than 1%) by the approximate formula 
of soft photon emission. 

A comparison of this plot with the corresponding Figure [3]for /i ^ e{'~f) 
conversion, which was obtained in the approximation of soft photon emission 
by Eq. ^ shows that the relative probabilities of these processes are quite 
close even though the energies released in these processes are somewhat 
different. 
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Conclusion 



This study has shown that radiative corrections including a virtual photon 
correction and soft photon emission below 1.5 MeV lead to a reduction by 
about 10% in the probability of fi ^ e conversion process calculated without 
radiative corrections. It is not necessary to consider a photon emission above 
1.5 MeV because these events do not survive a cut in the electron momentum 
of 103.5 MeV. 

The soft photons emission below 1.5 MeV contributes to a change of 
electron spectrum from monoenergetic one at 105 MeV to a spectrum with a 
low energy tail. The effect of smearing of the initial momentum distribution 
due to the soft photon emission is small in comparison with a smearing due 
to energy loss fluctuations in the target. The average energy of photons 
emitted below 1.5 MeV is found to be small of about 40 keV. 

By analyzing the radiative pion decay which possesses the same kine- 
matic properties and a comparable energy scale as ^ ^ 6(7) conversion it 
was found that up to the photon energy w = 5 MeV the exact photon spec- 
trum can be approximated by soft photon emission dN/duj ~ l/w with the 
precision better than 10%. 

We wish to thank A. Mincer and P.Nemethy for fruitful discussions and 
helpful remarks. 
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Appendix 

Spectra of the radiative pion decay 

According to [5j the double differential rate of pion decay including the 
inner bremsstrahlung with respect to photon energy u) and electron energy 
E is given by 

d'r _a r(vr ^ eu) ^^^^ 



dxdy 2ti 1 - A2 
where r(7r ^ ev) is the decay rate for non-radiative pion decay, A 
In Eq.([IO]) IiB is 



IlB 



1 - y + A2 



x'^{x + y — 1 — A-^) 



' 1 , .2. 2xA2(l-A2) 

x2 + 2(1-x)(1-A2)+ ^ ^ 



x + y - 1 - A2 

(11) 

where x = 2uj/m.,j and y = 2i?/m,r are relative photon and electron 
energies, respectively. These energies are directly measurable in radiative 
pion decay. 

Differential distributions in photon energy oj and electron energy E are 
obtained from Eq.f llOp by integrating over corresponding kinematic regions. 

It follows from the conservation of four-momenta p-,^ = k + q + pe, where 
PTj,k,q,pe are the four-momenta of pion, neutrino, photon, and electron, 
respectively, that 

k"^ = {j,^ - q -pg)2 = rn^ - 2mT,{ijj + E) + 2uj{E - \pe\z) + ml = 

where z = cosO , 9 is the angle between the electron and photon mo- 
menta. 

Since [z] < 1 the previous condition implies that 

{l + A'^-x-y + xyl2f - x'^{y'^/A - A^) < 0. 

The ranges of definition of x and y are easily established from this con- 
dition: 

2A < y < l + AM-y/2-(v/(y2-4A2))/2 <x< l-y/2+(v/(y2-4A2))/2 
or 
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< X < 1 - A^, (1 + - 2x + - x) < y < 1 + A^. 

Additional conditions were used to get ranges above: — 4A^ > 0, 
{Pn-Pef >0, (1 + A2-2x + x2)/(1-x) < 1 + A2. 

By integrating the double differential rate Eq.( fTOj) in corresponding 
limits one gets for the photon spectrum 



dr a r(7r ev) 



(x - 2(1 - A^ 



dx 2tt' 1 - A2 
) + 2(1-A2)Vx)M^) 



) -X(l- A2 



/(1-x)) 



and for the electron spectrum 



dV a r(7r — > ev) 1 



■{A + B) 



dy 2tt' 1 - A2 ■ 1 + A2 - y 



A = 2(1 - /\^){2/\^ -y)ln 



(l-y/2 + 7(y2_4^2)/2)2 

1 + A2 - y 



B = {1 + A2(-2 + 5A2) - 4A2y + y'^)ln 



(y/2-A2 + y(y2_4^2)/2)2 

A2(l + A2 -y) 
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